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Natural Quintessence and Large Extra Dimensions 
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We examine the late-time (nucleosynthesis and later) cosmological implications of brane-world 
scenarios having large (millimeter sized) extra dimensions. In particular, recent proposals for un- 
derstanding why the extra dimensions are so large in these models indicate that moduli like the 
radion appear (to four-dimensional observers) to be extremely light, with a mass of order 10~ 33 eV, 
allowing them to play the role of the light scalar of quintessence models. The radion-as-quintessence 
solves a long-standing problem since its small mass is technically natural, in that it is stable against 
radiative corrections. Its challenges are to explain why such a light particle has not been seen in 
precision tests of gravity, and why Newton's constant has not appreciably evolved since nucleosyn- 
thesis. We find the couplings suggested by stabilization models can provide explanations for both of 
these questions. We identify the features which must be required of any earlier epochs of cosmology 
in order for these explanations to hold. 
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I. INTRODUCTION 



If current experimental indications are to be believed, 
the universe is composed of islands of luminous material 
embedded amongst darker baryonic matter, wrapped in 
the riddle of nonbaryonic dark matter inside the enigma 
of dark energy. Best fits to cosmological models give the 
last three of these in the rough relative proportions of 
few% : 30% : 70% 0. It is staggering that essentially 
nothing is known about either of the two dominant com- 
ponents to the universe's present energy density. 

The very existence of dark energy presents serious puz- 
zles. The simplest account of the cosmic acceleration 
seems to be a non-zero cosmological constant A. However 
the energy scale of the cosmological constant is problem- 
atic, requiring an energy density pa ~ v with v ~ 10~ 3 
eV. The long-standing cosmological constant problem g] 
stems from the fact that vacuum fluctuations contribute 
to the cosmological constant, and simple estimates calcu- 
lated by assuming a Planck-mass cutoff give A a value 120 
orders of magnitude greater than is be consistent with 
observations. This discrepancy is reduced somewhat in 
models where the natural cutoff is the supcrsymmctry 
breaking scale, but even in that case there is a very seri- 
ous problem. 

Many physicists are attracted by the idea that the cos- 
mological constant problem could eventually be solved by 
some symmetry argument that would set A precisely to 
zero (although very interesting alternative views do exist 
P{. The observed cosmic acceleration creates problems 
with that line of reasoning, since it appears to indicate 
that A actually is non-zero. One way out is to assume 
that A really is zero, and that something else is causing 
the cosmic acceleration. This "something else" is usually 



called "dark energy" or "quintessence" 0|| . 

For quintessence, typically one proposes the existence 
of a scalar field whose mass is small enough for it to be 
evolving in a cosmologically interesting way at present. 
So far, essentially all quintessence models present the fol- 
lowing challenges or puzzles [pi. PI: Why don't radiative 
corrections destabilize the fantastic hierarchy between v, 
the scalar masses which are required {m ~ v 2 /M p ~ 
10 -33 eV) and the other known scales of physics? (Only 
the first quintessence model attempts to address this 
problem J4|) P2: Why isn't the very long-range force me- 
diated by the new light scalar detected in precision tests 
of gravity within the solar system? 

We here report on a class of quintessence models which 
is motivated by recent efforts to understand the hierar- 
chy problem, in which all observed particles except gravi- 
tons (and their supersymmetric partners) are trapped 
on a (3+l)-dimensional surface, or branc, in a higher- 
dimensional 'bulk' space J7],||. It turns out that within 
this framework the above two puzzles turn out to be re- 
solved in a natural way. We believe ours is the only 
quintessence models to successfully address both these 
particular challenges]]. 

A great virtue of the class of models which we explore is 
that they generally arise in the low-energy limit of theo- 
ries which were originally proposed to solve the gauge hi- 
erarchy problem, without any cosmological applications 



* Ref. M discusses another type of quantum correction which 
corrects classical rolling of the quintessence field off a local 
maximum. Those calculations do not apply to the type of 
potentials considered here. 



in mind. As we shall see, the generic low-energy limit of 
these models looks like scalar-tensor gravity [|o| , with the 
four-dimensional scalars of relevance to cosmology com- 
ing from some of the extra-dimensional modes of the six- 
dimensional metric. We show here that although these 
models were not proposed with cosmology in mind, the 
couplings they predict for these scalars have following 
properties: 

1. Good News: They can predict a late-time cosmol- 
ogy which provides a quintessence-style description 
of the dark energy. And they do so with scalars 
whose masses are extremely small, yet technically 
natural. 

2. Good News: They can evade the constraints on 
the existence of very long-ranged forces because the 
couplings of ordinary matter to the the very light 
scalars are predicted to be field dependent, and so 
evolve over cosmological timescales. As has been 
observed elsewhere pT|, once free to evolve, such 
couplings often like to evolve towards zero, imply- 
ing acceptably small deviations from General Rela- 
tivity within the solar system at the present epoch. 

3. Good News: The expectation values of the scalar 
fields whose cosmology we explore determine the 
value of Newton's constant, and the success of Big- 
Bang nucleosynthesis (BBN) implies strong lim- 
its on how much this can have changed between 
then and now. We find the evolution of New- 
ton's constant changes very little over cosmologi- 
cally long times, for a reasonably wide range of ini- 
tial conditions when we enter the BBN epoch: New- 
ton's constant tends not to change once the motion 
of the scalar field is either kinetic- or potential- 
dominated. 

4. Bad News: The previous two cosmological suc- 
cesses are not completely generic, in the following 
sense. First, they depend on dimensionless cou- 
plings taking values which are of order e ~ 1/50 or 
so. (Although couplings of this order of magnitude 
naturally do arise in the models of interest, we re- 
quire some modest coincidences in their values to 
ensure sufficiently small evolution in Newton's con- 
stant.) Second, although many initial conditions 
do produce an acceptable cosmology, this success 
is not completely generic. Attractor solutions ex- 
ist which would cause too much variation in New- 
ton's constant if they described the cosmology be- 
tween BBN and today. Furthermore, obtaining the 
present-day value of Newton's constant requires ad- 
justing the features of the scalar potential. 

We present our results here as an example of a vi- 
able cosmology, based on using the brane-world's extra- 
dimensional 'radion' mode as a quintessence field. Our 



ability to do so may itself be considered something of a 
surprise, since extra-dimensional aficionados had earlier 
examined the radion as a quintessence candidate, and 
found it wanting. Our conclusion here differs, because of 
our recognition that the radion couplings tend to be field 
dependent, in a way we make precise below. 

We regard this work as a first exploration of the con- 
ditions which a successful cosmology must assume. We 
believe our results to be sufficiently promising to justify 
further exploration of this class of models, to see to what 
extent the unattractive features can be improved upon, 
and to ascertain what observations can test whether the 
universe indeed proceeds as we describe here. 

Our presentation is organized as follows. In the next 
section we briefly describe the brane models we use, and 
their motivations and problems. The low-energy action 
relevant for cosmology is also presented in this section. 
Section 3 describes some of the cosmological solutions of 
the theory, focusing in particular on analytic expressions 
on which intuition can be based. More complicated, and 
realistic, solutions are then described numerically. We 
close in Section 4 with a discussion of our results. 



II. THE MODEL 

Our point of departure is the brane-world scenario 
in which all observed particles except the graviton are 
trapped on a (3+l)-dimensional surface, or brane, in a 
higher-dimensional 'bulk' space. In particular, we imag- 
ine the scale of physics on the brane is Mb ~ 1 TeV and 
the bulk space is approximately six-dimensional, com- 
prising the usual four, plus two which are compact but 
with radius as large as 1/r ~ 10~ 3 eV [0,0 • The large ra- 
tio Myr ~ 10 15 is phenomenologically required to repro- 
duce an acceptably weak gravitational coupling in four 
dimensions, M p = (SttG)- 1 / 2 - M fc 2 r - 10 18 GeV. 

These kinds of models face two main difficulties, one 
theoretical and one phenomenological. Theoretically, 
these models trade the large hierarchy Mb/M p — 10~ 15 
for the hierarchy Mj,r = 10 15 , but an explanation of 
this last value ultimately must be provided. Recently 
progress towards providing this explanation has been 
made, by proposing energetic reasons why the extra- 
dimensional radii might want to prefer to take large val- 
ues p2lfl| ■ O ur quintessence model does not rely on the 
details of these mechanisms, but our choice for the scalar 
field couplings is strongly motivated by the mechanism 
ofref. g| . 

Phenomenologically, there are strong upper limits on 
how large M\,r can be. These limits come from super- 
nova Mjl5| and cosmological observations Cm which 
preclude the existence of Kaluza-Klein (KK) partners 
of the bulk-space graviton and its superpartners. We 
believe the present-day observational constraints on KK 
modes can be accommodated if only two dimensions have 



r as large as (10 _3 eV) _1 , but prevent r from being taken 
any larger. Although the cosmological constraints gener- 
ally impose stronger upper limits on r they depend more 
on the details of the assumed cosmology, and can be 
evaded under certain conditions (which we describe in 
more detail later). 



A. The Low-Energy Action 

We start by sketching the brane world scenario, and 
its low energy action. Our discussion closely follows that 
of ref. |J|. 

The couplings of the brane-world model relevant to 
cosmology are those of the six-dimensional metric, Gmn- 
(Six dimensions are required in order to allow the ex- 
tra dimensions to have radii as large as we need.) The 
leading terms in the derivative expansion for the six- 
dimensional action have the form S = Sb + £fc, where 
the Bulk action describing the metric degrees of freedom 
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where 1Z denotes the scalar curvature built from the six- 
dimensional metric, and we write no cosmological con- 
stant for the bulkJT] 

The dependence of the brane action on its various mat- 
ter fields is not crucial for our purposes, but what is im- 
portant is that it depends only on the metric evaluated 
at the brane's position: 



S b = Id Xy/gC b (g 
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where the ellipses denote all of the other fields on which 
Cb could depend. 

The tree-level dimensional reduction of this action us- 
ing the metric 



fMN 



9nv(x) 

p 2 {x)h mn (y) J ' 



(3) 



with p = Mb r, gives the leading contribution to the ef- 
fective four-dimensional Lagrangian: 
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^This can be ensured by symmetries, such as in a supersym- 
metric compactification fl?]]. Like everyone else we assume 
the vanishing of the effective four-dimensional cosmological 
constant, and our ideas do not add to the understanding of 
why this should be so. 



where we have adopted the conventional normalization 
Ml J d 2 y \fh = 1. Here R(g) and R(h) denote the cur- 
vature scalars computed from the metrics g^ v and h mn , 
respectively. 

Quantum corrections modify this action in several im- 
portant ways. First, the kinetic terms can become renor- 
malized, introducing the possibility that their coefficients 
can acquire a dependence on r. Second, terms not writ- 
ten in eq. (||) can be generated, including in particular a 
no-derivative potential, U(r). At large r, U(r) can typ- 
ically be expanded in powers of 1/r, with the leading 
term falling off like l/r p , for some p > 0. For instance, 
for toroidal compactifications, the leading contribution 
can arise due to the Casimir energy, which produces the 
potential U(r) = Uo/r 4 H . 

We are in this way led to consider a scalar-tensor 
theory [ |10|Jll]j20(| , for the couplings of r and the four- 
dimensional metric, which we write as follows: 
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As just described, the tree-level dimensional reduction 
of the 6d Einstein-Hubert action predicts A = B = 1 and 
U = 0, and with these choices eq. (||) describes a Jordan- 
Brans-Dicke theory with coupling parameter lo = —1/2 
fll8| . Notice that the matter Lagrangian, C m , describes 
the couplings of all of the known particles, and only de- 
pends on g^y but not on r. (Ordinary matter does not 
couple directly to r in the brane-world picture because it 
is trapped on the brane, and does not carry stress energy 
in the extra dimensions.) 

Notice that if the scale of U(r) is completely set by r, as 
is true for the cases discussed in ]13| , then the mass which 
is predicted for this scalar is incredibly small: of order 
to ~ v 2 /M p ~ 10 -33 eV. It is the existence of such a small 
mass which allows r to potentially play a cosmologically 
interesting role right up to the present epoch, and so to 
be a candidate for a quintessence field. 

Since the predictions for A, B and U depend crucially 
on the approximation made in their derivation, in general 
we expect them to be given by more complicated func- 
tions. This turns out to be true, in particular, within the 
hierarchy-generating mechanism described in ref. p3[ . 
Since we use this proposal to motivate our choices for 
these functions, we briefly pause here to summarize its 
implications. 

The proposal supposes that the spectrum of bulk 
states includes a six-dimensional scalar, 0, which counts 
amongst its interactions a renormalizable cubic self- 
coupling, g<f) 3 . Since the coupling g of this interac- 
tion is dimensionless, it can introduce corrections to the 
low-energy action which are logarithmic in r. Since all 
other couplings have negative mass dimension in 6d, their 



quantum contributions are suppressed by powers of 1/r, 
implying the coupling g can dominate at large r. 

The bottom line of this scenario is that the domi- 
nant corrections to A, B and U are of the form 1 + 
c\ a(r) + C2 a 2 (r) + ■ ■ ■, where the running loop-counting 
parameter is given by a(r) = ao + aQlog(r/ro) + • • ■, 
with «o = 9o/(4 7r ) 3 - I n this scenario the potential 
U(r) is plausibly minimized for a(r) ~ 1, implying a 
minimum for which there is the desired large hierarchy: 
Mi>r ~ exp[+l/ao] 3> 1. It is generic in this scenario that 
the quantum corrections to the other low-energy func- 
tions, A and_B, are also likely to be large. 

For illustrative purposes when exploring cosmology, 
we choose the functions A, B and U which are sug- 
gested by this picture in the perturbative regime, for 
which A « 1 + alog(M b r), B w 1 + blog(M b r) and 
U sw (^r) [1 + clog(M&r)]. Here a,b and c are small, 
since they are proportional to the dimensionless coupling 
constant, e ~ g 2 /(Air) 3 , and higher orders in e give higher 
powers of log(M b r) although these have not been writ- 
ten explicitly. In the spirit of this picture we suppose 
e to be of order 1/50, as required in order to ensure 
that U is minimized when M^r ~ 10 15 . For cosmological 
purposes, we shall be interested in the regime for which 
elog(M;,r) < 1 is negligible, but elog 2 (Af b r) is not. We 
take Uq to be generic in size: Uq 
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B. Naturalness 

Since naturalness is one of our motivations, we briefly 
summarize the arguments presented in ref. |l3fl as to why 
radiative corrections do not destabilize the scalar mass. 
The argument takes a different form for the contribution 
due to the integrating out of scales above, or below, v ~ 
1/r ~ 1CT 3 eV. 

For scales below v, the effective theory consists of a 
four-dimensional scalar-tensor theory coupled to itself, 
and to ordinary matter with couplings, K, of gravita- 
tional strength: k ~ l/M p . Radiative corrections to the 
scalar mass within this effective theory are then generi- 
cally acceptably small, being of order 6m 2 ~ k 2 A 4 , with 
ultraviolet cutoff A ~ v. These are therefore the same 
size as the tree-level mass, as is required if they are to be 
technically natural. 

Radiative corrections to the masses coming from higher 
scales, between v and M&, are the dangerous ones. Gener- 
ically, these are too large and it is the integrating out 
of these scales which introduces the naturalness prob- 
lem to most quintessence models. Within the framework 
of ref. |l3[ , the contributions from these scales are sup- 
pressed by supersymmetry, which is proposed to be bro- 
ken at the scale M b on the brane on which observable 
particles live. (This ensures acceptably large splittings 
between the masses of ordinary particles and their su- 
perpartners.) 



Supersymmetry suppresses corrections to the radion 
mass because it relates r to the four dimensional gravi- 
ton, which is massless. Indeed, above the scale v the bulk 
space is really six-dimensional, and the radion is simply 
a component of the six-dimensional metric. As such, its 
properties are tied to those of the four-dimensional metric 
by the supersymmetry of the bulk space. Since the metric 
and its superpartners only couple to the supersymmetry 
breaking with gravitational strength, supersymmetry- 
breaking effects in the bulk sector are suppressed to be 
of order M%/M p ~ 1/r ~ v. 

Since the scalar potential for the radion cannot arise 
until supersymmetry breaks, it is naturally of order 
v 4 F(ip) where ip is the dimensionless component of the 
six-dimensional metric which defines the radion. tp'a 
kinetic term, on the other hand, comes from the 6d 
Einstein-Hubert action, and so is allowed by unbroken 
supersymmetry. At tree level it is of order M 2 (d(p) 2 . 
and the coefficient M 2 is much larger than the contri- 
butions produced by integrating out physics at scales 
smaller than Mb <C M p . Again we are led on dimen- 
sional grounds to a mass which is of order m ~ v 2 jM p . 



C. Cosmological Equations 

The cosmology of the radion is most easily described by 
passing to the Einstein frame, for which the radion and 
metric are canonically normalized. For the choices of A = 
1 + alog(Mfer) and B = 1 + 61og(M b r) the canonically 
normalized variables are: 



g^=g^/[{M b r) 2 A 2 



and log(Affcr) = a%(l — b\), 

(6) 



with a = (l-3o/4)/2, b = (a-6)/16 and \ = x/M b - This 
choice of variables is called the Einstein frame, to distin- 
guish it from the Jordan frame, as defined by eq. (g)Jj 
In the Einstein frame the action, eq. (|5j), becomes 
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exp[-A x(l -bx)], 



'A note on numerics: In the Jordan frame M p ~ M^r, ordi- 
nary particle masses are m p ~ Mb and Kaluza-Klein masses 
are m^k ~ 1/r, while in the Einstein frame M v ~ Mb, 
m p ~ 1/r and m^ ~ l/(M(,r 2 ). Clearly, ratios of masses 
like mkk/m-p = l/(Mf,r) are the same in either frame, but 
Mb — 1 corresponds to using TeV units in the Jordan frame, 
but Planck units in the Einstein frame. 



and Ao = 8a. (Even though b = 0(e) we do not expand 
the second term in the exponent in this expression be- 
cause such an expansion presupposes ex 2 <§C 1, which is 
not true for our applications.) 

For cosmological applications we take ds 2 = — dt 2 + 
a 2 dx 2 in the Einstein frame, in which case the Lagrangian 
(pi) implies the usual Friedman equations: 
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(9) 



where both the energy density, p, and pressure, p, receive 
contributions from x and the radiation and matter sec- 
tors - p = p r + p m + p x and p = p r + p m + p x , where 

p r ~ \Pr and p m f=S 0. 

In the Einstein frame the matter contributions also 
depend explicitly on %: 
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where the quantities p m o and p r Q, may be related to 
the density of matter and radiation at the present epoch 
(which we choose to be a = 1). 

Eqs. (||) must be supplemented with the evolution 
equation for %, which is: 



X + 3Hx + V'(x) + 
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where T^ v is the stress energy of the radiation and mat- 
ter, and T](x) = — o(l + a/2 — 2o%). The successful com- 
parison of General Relativity with precision tests of grav- 
ity within the solar system require ?7 2 (xo) Ss 10~ 3 when 
evaluated at the present epoch, x — Xo [0pQ]- Notice 



that since radiation satisfies T 7 ^ = 0, only matter con- 
tributes to this last coupling, which is proportional to 
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Pn 



III. COSMOLOGY 

We now show that cosmologically reasonable solutions 
may be found to eqs. (||) and (|ll|), which satisfy all of 
the bounds to which the theory must be subject. Since 
our goal is to demonstrate the existence and illustrate 
some of the features of these solutions, our analysis is not 
systematic and should be regarded as only preliminary 
explorations of models of this type. 



A. Constraints 

Our main assumption is that the universe enters into 
an effectively four-dimensional radiation-dominated evo- 
lution at some temperature T» higher than a few MeV, 
and so before Big Bang Nucleosynthesis (BBN) . This as- 
sumption is a prerequisite for describing cosmology in 



terms of four-dimensional fields, and seems likely to be 
required if the success of standard BBN is not to be for- 
saken. 

There is a limit to how high T* may be chosen, since the 
four-dimensional picture requires the extra-dimensional 
KK modes to be cold. If T* is too high, then these KK 
modes can be thermally excited too efficiently through 
their couplings with ordinary particles. If these modes 
are too abundant they can overdose the universe or decay 
too frequently into photons after recombination, and so 
contribute unacceptably to the diffuse gamma ray back- 
ground J7||l6|]. If the KK modes all decay into observ- 
able daughters, then these constraints require T» not to 
be larger than the TeV scale, and can preclude having 
1/r ~ 1CT 3 eV. 

As the authors of refs. [f7|Jl6[ themselves point out, 
these constraints may be evaded, depending on the cos- 
mology which is assumed at earlier times in the universe's 
history. For instance, contributions to the diffuse gamma 
ray background is suppressed if the KK modes decay 
quickly enough into unobserved modes to be no longer 
present after the recombination epoch, when photons de- 
couple from matter. Such fast decays can be arranged, 
for instance, if the KK modes decay more frequently into 
particles on other branes. Such decays may, in fact, be 
likely in cosmological scenarios for which the very early 
universe is dominated by a gas of branes and antibranes 
pl[ . (Another escape is possible if the extra-dimensional 
space is hyperbolic instead of toroidal, with its volume 
much greater than the appropriate power of its radius of 
curvature: V ^> r 2 |22]. In this case KK modes can be 
very heavy, m^k ~ 1/^c, even though the hierarchy prob- 
lem is solved by having V large.) Although our present 
interest is for T* < 1 TeV, we have no difficulty imagining 
some variation on these themes being invoked if scenarios 
requiring T* > 1 TeV were of interest. 

Our second major constraint is the requirement that 
the quintessence field, Xi n °t destroy the successes 
of BBN. This constraint arises because the success of 
BBN may be regarded as indicating that the hierarchy 
M w /M p ex (Gn/Gf)^ 1 ^ 2 is not much different at the 
epoch of nucleosynthesis than it is today. A sufficient 
condition for this to be true is to ensure that M b r not to 
have evolved by more than about 10% of its present value 
between BBN and the present epoch. This constraint is 
much stronger than what is required by garden-variety 
quintessence models. It is very strong because it requires 
that the field r cannot have rolled significantly over cos- 
mologically long times. A special property of the cosmo- 
logical solutions we shall examine is that they often have 
the property that r becomes fixed over cosmologically 
large timescales. Motivated by the BBN constraints we 
are particularly interested in those cosmological solutions 
for which r does not strongly vary after BBN. We shall 
see that such solutions can arise for reasonably generic 
initial conditions for the functions V(%) and T}(x) we are 



using. 

Next, we must demand agreement with present-day ob- 
servations, which come in several important types. First, 
precision measurements of gravitational forces indicate 
that the coupling r](x) must be < 0.03 in order not to be 
in conflict with the success of General Relativity. Sec- 
ond, if x is to provide an explanation of the dark energy, 
then its energy density must now be beginning to domi- 
nate that of matter, with w = P x l Px < ~ I ^ this ener gy 
density is to be currently accelerating the universe. 

Finally, we require \ n °t to ruin the success of scale- 
invariant fluctuations in describing present-day obser- 
vations of the cosmic microwave background radiation 
(CMB). Since a detailed analysis of the CMB in these 
models is beyond the scope of the present paper, we in- 
stead apply a poor-man's bound by requiring cosmolog- 
ical properties not to be too different from the standard 
picture near the epochs of radiation/matter crossover and 
recombination. 



B. Approximate Cosmological Solutions 

Before laying out a cosmology which satisfies the above 
constraints, it is worth building intuition by outlining 
several approximate analytical solutions to the cosmo- 
logical equations which turn out to control the features 
which are seen in the numerical results. Since the micro- 
physics motivates using the potential: V(x) — V p e~ Xx , 
with V p (x), -Mx) an d the coupling r)(x) varying much 
more slowly as (polynomial) functions of x than does the 
exponential e Xx , we explore solutions to the equations 
when A, V p and 77 are all constants. The cosmological 
evolution of such fields in the limit 77—5-0 has been much 
studied |23(| , and our results go over to these in the ap- 
propriate limits. 

There are two types of solutions in this limit which are 
of principal interest. These solutions may be classified 
according to which the terms dominate in the x equation, 
eq. (|ll). 



1. Scaling Solutions 

The first class of approximate solutions we consider 
are the scaling, or tracker, solutions, of which there are 
two different types. These solutions are characterized by 
having all energy densities scaling as a power of the scale 
factor, and either: (i) x ~ 3i?x ~ V ^ VPm or, (ii) 
X ~ 3H x ~ VPm ^> V . Both types of solutions can arise 
regardless of whether it is the scalar, radiation or matter 
which is dominating the evolution of the universe. 

Assume, therefore, the scale factor scales with time 
as a/ao = (t/to) k , and so H — a/a — k/t. From the 
Friedman equation, H 2 oc p, we see that if the dominant 



energy density scales as p oc a ™, then we must have 
k = 2/n. 

(i) x ~ 3-ffx ~F'» rjp m : In this regime we assume 



V oc e~ A * 



(12) 



and we find from the x field equation, eq. (11), that s = 



If radiation or matter dominate the energy density, 
then n = 3 or n = 4. In this case we find solutions 
for which K — ^x 2 oc V(x) oc 1/a™, and so r oc a™/ 8 . 
Eq. (O) also implies that these scaling solutions are 
trackers in the sense that p x tracks at precisely the frac- 
tion n/X 2 of the dominant energy density as the uni- 
verse expands. Solutions of this type only are possible 
if A 2 > n. For X 2 < n the scalar energy density tends 
towards scalar domination of the universal expansion. 

If it is x itself which dominates the energy density, then 
k = 2/X 2 and a scaling solution of the type we seek exists 
so long as X 2 < 6. In this case both K and V scale as 
\/a x , and so r oc a A / 8 . 

These are the tracker solutions which are familiar from 
standard discussions of scalars evolving with exponential 
potentials |23|,|24|. Clearly these tracker solutions imply 
a relatively rapid variation of r as the universe evolves. 

(ii) x ~ 3-ffx ~ T)p m ^> V: We next consider the case 
where the rfT^ M term dominates the scalar potential term 
in eq. (|ll|). Notice that this can hold even if both p r 
and p x are much larger than p mi so the tracker solution 
we here obtain can be relevant even when nonrelativistic 
matter makes up only a very small part of the universe's 
total energy budget. 

Under these assumptions, because the last term of 
eq. ( |ll| ) dominates the potential term, we ask p m to scale 
like x an d Hx'- 



Pm oc -s- oc — 
a°r a* 



(13) 



which, when inserted into eq. (|1L) implies s = n and so 
K oc p m oc p oc l/a n . As usual, if the universe is radia- 
tion or matter dominated then n = 4 or n = 3. (Inter- 
estingly, both matter and radiation fall as 1/a 4 within 
this kind of tracker solution in a radiation-dominated 
universe.) Alternatively, K oc p m oc p oc l/a A , if the 
universe is ^-dominated. 

Although the scaling of scalar kinetic energy is the 
same for tracker solutions of type (i) and type (ii), they 
differ markedly in the scaling which they imply for r. In 
the present solution the scaling of r vs a follows from 
the known scaling of p m . For radiation/matter dom- 
ination, where p oc 1/a™ we find ra 3 oc a™ and so 
r oc a"" 3 . Clearly r oc a for radiation domination, but 
r only varies logarithmically during matter domination! 
During a scalar-dominated phase, on the other hand, we 
have r oc a A _3 . Again, r does not vary if A 2 = 3. 



What is striking about these solutions where r changes 
so slowly is that this happens in spite of its not sitting at 
the minimum of the potential. It happens, even though 
the radion's kinetic energy dominates its potential en- 
ergy, because the amount of roll is limited by the Hubble 
damping generated by the dominant energy component 
of the universe. 



2. Transient Solutions 

What is perhaps even more striking than the trackers 
for which r slowly varies is the behavior of r during the 
transient evolution as the tracker - which is an attractor - 
is approached. Since the approach to the tracker solution 
may be approached from initial conditions for which K S> 
V or K -c V, we consider each of these in turn. 

K < V (Slow Roll): If K < V then it is a good ap- 
proximation to neglect \ in eq. (|llj) relative to the other 
terms. If r/p m is also negligible relative to V' , then the 
evolution is described by 3Hx + V = 0, which has as 
solutions: 



V(a) 



2V; x ita^ xit 



aP 



a? ■ 

"exit 



(14) 



where p = n for radiation (n = 4) and matter (n = 3) 
domination, and p = A 2 for scalar domination. Along this 
trajectory the kinetic energy grows as K/V = X 2 V/(6p). 
This solution implies V « V cx n is a constant (and so 
also must be r) until a ~ a C xitj after which point the 
solution joins the tracker solution, for which p x ex 1/aP 
and p x /p = p/\ 2 . 

K 3> V (Kinetic Domination): UK dominates the scalar 
evolution, then V' and r\p m may be dropped relative to % 
and 3Hx- In this case the x equation may again be solved 
analytically, with solution: x — Ci ~ C2 (do/ a) , 

where C\ , Ci are integration constants with C^ related to 
the initial velocity by C2 — (3 — l/p)ao(dx/da)o. Here, 
as usual, p = 2/n for radiation/matter domination, and 
p = 2/A 2 for scalar domination. 

This solution implies the kinetic energy, K = i% 2 
varies as: 



K(a)=K (^)\ 



(15) 



regardless of the value of p. 

The evolution of the potential energy is more interest- 
ing in these solutions. Consider first the case of radia- 
tion/matter domination, for which p = 2/n so 



V(a) = Vq exp < c 



ao\ 3 ~"/ 2 



(7) 



(16) 



As the universe expands according to this solution, x 
carries an ever-smaller fraction of the energy density since 
K(a) ex 1/a 6 falls faster than does the dominant matter 
or radiation, p ex 1/a™. Eventually K becomes compa- 
rable in size to V(a), which for large a asymptotes to a 
constant value: V(a) — + Voo — Voe~ c . At this point the 
approximation of dropping the scalar potential in eq. (O) 
breaks down, and we match onto the slow-roll solution 
discussed above. The amount of expansion which oc- 
curs before V and K are comparable is determined by 
a/a = (e c K /V )^ e . 

What is remarkable about this solution is that the con- 
stant c cannot be made arbitrarily large. To see this, 
notice that c can be related to the initial kinetic energy 
by using K = ^H 2 a 2 (dx/da) 2 . Using the upper limit 
K < 3-ff 2 , which is dictated by the inequality K < p and 
the Friedman equation, we see that c < \/6A(3 — 1/p). 
For A = 4 and p = 1/2 (radiation domination) we have 
c<4\/6~9.6. 

This upper limit on c implies V cannot fall by more 
than a factor of e~ c < e~ A ^ 6 during a radiation domi- 
nated epoch, despite the fact that r's evolution is being 
dominated by its kinetic energy. Since V ex 1/r 8 we see 
that r/ro < e A ^ 6 / 8 = e°' 3A varies hardly at all during 
the entire kinetic roll. Physically the paradoxical upper 
limit on r during a kinetic dominated roll arises due to 
the Hubble damping which increases with the kinetic en- 
ergy, if the kinetic energy gets large. This observation 
has also been made in another context in ref. |2q| . 

A similar result holds for the kinetic roll in the scalar- 
dominated phase, although with results which differ 
strongly depending on the value of A. If A 2 < 6, the 
above analysis applies directly, with p = 2/A 2 . As before 
K ex 1/a 6 and 



V{a) = Vq exp < c 



er /2 - 



(17) 



with c = (3 — l/p)Aao(<ix/^ a )o- 



Provided A 2 < 6, V(a) again asymptotes to Voe~ c , with c 
once more subject to a very general upper bound. Once 
K falls to be of order V, this solution matches onto a 
slow roll solution, and approaches the tracker solution, 
which exists for A 2 < 6. 

If A 2 > 6, however, the large-a limit of V(a) is very 
different. In this case V(a) ex exp[— \c\(a/ao)( x ~ 6 - )/ ' 2 ], 
implying V(a) falls very quickly to zero as the roll pro- 
ceeds. This solution describes the attractor for the case 
A 2 > 6, which is a kinetic-dominated roll rather than a 
tracker solution. 



3. Realistic Complications 

Although these solutions describe well most features 
of the numerical evolution we describe below, there are 



also important deviations. These deviations arise be- 
cause rj and the quantities A and V p (of the potential 
V = V p e~ Ax ) used in the numerical evolution are all 
polynomials in \- Motivated by the microphysical dis- 
cussion given in section II we choose: 



v = m + mx; Hx) = A + Aix; 



(18) 



existence proof that cosmologically interesting evolution 
can really be found using the scalar properties which are 
suggested by the microphysics. We find that the main 
constraint which governs the construction of such solu- 
tions is the requirement that r not vary appreciably be- 
tween BBN and the present epoch. 



and 



V P = Vo 



V 2 O 



(19) 



with tjq, X and Vo all taken to be 0(1); ?7i,Ai and V\ 
assumed to be 0(e) and V 2 taken 0(e 2 ), where e ~ 1/50. 
(The second-order quantity V 2 is included here in order 
to ensure V does not change sign for any value of %.) 

Because these quantities vary more slowly with \ than 
does the exponential in V, their evolution over cosmolog- 
ical timescales is also slow, and the above approximate 
solutions provide a good description for most of the uni- 
verse's history. We have found that their x-dependence 
can nevertheless play an important role when construct- 
ing viable cosmologies, by allowing important deviations 
from the approximate solutions, including: 

• The %-dependence of V p becomes important for 
those x satisfying Vp/V p ~ A, since the potential 
V can acquire stationary points near such points. 
These kinds of stationary points can slow or trap 
the cosmological \ rou i m much the same way as 
was exploited in ref. p6| . 

• The %-dependence of A can cause transitions be- 
tween the domains of attraction of the various at- 
tractor solutions. For instance, in a scalar dom- 
inated era, the evolution of \ can drive A across 
the dividing point, A 2 = 6, forcing a transition be- 
tween the scaling and kinetic-dominated attractor 
solutions. 

• The %-dependence of rj is crucial for evading the 
bounds on long-range forces during the present 
epoch, because it allows r\ to vanish for some choices 
of x- The bounds can be evaded provided the cur- 
rent value of x IS sufficiently close to such a zero of 
77. Furthermore, as described in detail in ref. pi| , 
in many circumstances x' s cosmological evolution 
is naturally attracted to the zeros of 77 as the uni- 
verse expands. 

Because of the natural hierarchy satisfied by the co- 
efficients of these functions, the consequences of x~ 
dependence generically occur for x ~ 0{l/e), corre- 
sponding to r close to its present value. 



C. A Realistic Cosmological Evolution 

We next describe a realistic cosmology which exploits 
some of these properties. We present this model as an 
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FIG. 1. The logarithm (base 10) of the energy density (in 
Planck units) of various components of the universe, plotted 
against the logarithm (base 10) of the Einstein frame scale 
factor, a (with a = 1 at present). The various curves show 
the energy density of matter (dashed), radiation (dotted) and 
scalar energy (solid). The dot-dashed curve shows the scalar 
kinetic energy density, which coincides with the total scalar 
density except at late times 

In order to fix ideas, we imagine starting the universe 
off at T„ > lMeV, perhaps after an epoch of earlier 
inflation, in a kinetic-dominated roll. We note that such a 
state is the generic endpoint if we start with A ~ 4 > \/6, 
as is motivated by the tree-level expression for the radion 
potential, and assume the scalar dominates the energy 
density of the universe. (As we will discuss briefly later, 
one might chose to make other assumptions about the 
initial conditions.) 

With this assumption the scalar energy density must 
eventually fall below the density of radiation, thus initi- 
ating the radiation-dominated era. We assume that this 
radiation-domination begins shortly before BBN, so that 
nucleosynthesis occurs during radiation domination. 

Our second assumption is dictated by the requirement 
that r does not change appreciably between BBN and 
now. Although the kinetic-dominated roll before radia- 
tion domination will automatically ensure that F«i(, 
we ask that V be within an order of magnitude of to its 



present value as radiation domination begins. This en- 
sures that r will have an initial value at the onset of ra- 
diation domination which is very close to its present-day 
value. We make this (fairly ugly) assumption to demon- 
strate the feasibility of models along the lines we are pre- 
senting, but hope to be able to eventually identify a more 
natural frameworks in which such a condition could be 
ensured by the system's dynamics rather than as an ini- 
tial condition. 
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FIG. 2. The same plot as Figure QJ, but with the energy 
densities given as fractions of the critical density. The various 
curves show the energy density of matter fl m (dashed), radi- 
ation Q r (dotted) and total scalar energy f2 x (solid), while 
the dash-dot curve gives the equation of state parameter, 
w = p/p, for the scalar field. 

For A ~ 4, once the universe becomes radiation domi- 
nated the scalar in the domain of validity of the tracker 
solutions described above. With the given initial condi- 
tions the scalar enters radiation domination in a kinetic- 
dominated roll, and so tries to approach a tracker solu- 
tion by first damping its kinetic energy in a continued 
kinetic-dominated roll, and then performing a slow-roll 
at fixed r. We have seen that the universe expands by 
an amount (a/oo) 2 ~ Kq/Vq during the transient phase, 
before the tracker solution is reached, so given the low 
values of Vq assumed at the onset of radiation domina- 
tion, Kq/Vq is very large, and so the transient evolution 
to the tracker can easily take the entire time from BBN 
to the present epoch. 

We have seen on general grounds that the value of 
r does not change by more than a factor of 0(1) dur- 
ing this transient phase, including both the transient's 
kinetic-dominated and slow-rolling parts. The challenge 



in making this scenario realistic is to have this transient 
evolution survive right up to the present day, without 
first being intercepted by a tracker solution. 
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FIG. 3. The logarithm (base 10) of the dimensionless quan- 
tity Mj,r, plotted against the logarithm (base 10) of the Ein- 
stein frame scale factor, a (with a = 1 at present), showing 
how the radion has not evolved appreciably between nucle- 
osynthesis and the present epoch. 

We have been able to find initial conditions and plau- 
sible couplings for which this is accomplished. One such 
is illustrated in Figures (1) through (4), which plot the 
evolution of various quantities from the onset of radiation 
domination to the present. The evolution described by 
this simulation uses the values suggested by the lowest- 
order microscopic action, 770 = 1, Ao = 4 and Vq = 1, 
plus the following correction terms: 771 = —0.015, Ai = 0, 
Vi = -0.030 and V 2 = 0.00046. 

Figure (Tfl) shows how the energy density in the uni- 
verse is distributed between radiation, matter and the 
kinetic and potential energy of the scalar field, as a 
function of the universal scale factor, a, in the Einstein 
frame. Figure (0) plots the same information, but nor- 
malized as a fraction of the total universal energy density, 
and also plotting the scalar equation-of-state parameter, 
w = p x /p x - Figure (||) gives a plot of r vs scale factor, 
showing that r does not vary appreciably between BBN 
and the present epoch. Finally, Figure (|J) shows the evo- 
lution of 77 against scale factor, showing that 77 evolves to 
sufficiently small values at the present epoch. 

As is evident from the figures, this solution satisfies 
all of the cosmological and present-day bounds we were 
listed earlier. This was accomplished by choosing appro- 
priately the constants governing the scalar interactions, 



by ensuring V to have a minimum relatively near to a zero 
of 77. This was required because although r naturally does 
not evolve significantly for most of the universe's history 
(after the onset of radiation domination), for pure expo- 
nential potentials it generically tends to join a tracker 
solution near the onset of matter domination. Because r 
evolves too far once in this tracker solution, it must be 
avoided. The solution illustrated in the figure does so by 
having \ become snagged by the minimum of V before 
entering into the matter-dominated phase. 
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FIG. 4. The scalar coupling function, rj, plotted against the 
logarithm (base 10) of the Einstein frame scale factor, a (with 
a — 1 at present), showing how r\ is small enough to satisfy 
the present-day bounds on the existence of long-range scalar 
forces. 



IV. DISCUSSION 

We have shown that viable cosmologies may be built 
using the radion as the quintessence field. We use ra- 
dion interactions which were suggested by a microphysi- 
cal model which was recently proposed to naturally solve 
the hierarchy problem within a brane-world framework. 
The radion in this model is naturally light enough to 
play a role in late-time cosmology, and this small mass 
is technically natural in that it is not destabilized by 
radiative corrections. The model also satisfies all cur- 
rent bounds which constrain the existence of long-range 
scalar-mediated forces. 

The success of the radion in filling the quintessence role 
is somewhat of a surprise, since radion models have long 
been believed to have phenomenologically unacceptable 
couplings. The model evades these problems because the 



couplings are predicted to be weakly field dependent, al- 
lowing the scalar couplings to evolve over cosmological 
timescales. The model evades all present-day bounds on 
new forces because the relevant couplings evolve towards 
small values at late times. 

Although most of the parameters and initial conditions 
chosen were natural in size, the model does have two fea- 
tures which we believe need further improvement. First, 
although all scalar couplings were chosen with natural 
sizes, the precise values chosen were adjusted to arrange 
the scalar potential to have a minimum close to a zero 
of r\. We would prefer to find an attractor solution, a la 
Damour and Nordtvedt jll| , which more naturally draws 
77 to small values in the present epoch. 

A second unsatisfactory feature of the cosmology pre- 
sented was its reliance on the entering of radiation domi- 
nation with r having close to its present value. Although 
large values of r are naturally generated by kinetic- 
dominated rolls during scalar dominated epochs, there 
is no natural reason why they should have precisely the 
current value. This problem might be resolved by a 
more developed understanding of radion initial condi- 
tions. Perhaps one could argue that there are no a-priori 
constraints on the initial conditions for r. In that picture, 
it may be the case that the largest amount of phase space 
that matches the observation that Newton's constant is 
constant and non-zero today would be one where r has 
the current value throughout the cosmic evolution. (In 
that reasoning, the constancy of G is not a prediction, 
but an input.) 

A natural way to alleviate both of these problems 
would be to consider more moduli than just the ra- 
dion, since other moduli would be less constrained by 
the requirement that Newton's constant not appreciably 
change since nucleosynthesis. We believe work along this 
lines to be worthwhile, given the model's successful ad- 
dressing of the naturalness issues of quintessence models, 
and the promising cosmology to which it leads despite 
its being devised to solve purely microphysical problems 
such as the hierarchy problem. 
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